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Within the semiclassical approach, we calculate contributions of the Berry phase and of the 
Zeeman coupling of an electron moment with the magnetic field to the phase of the Shubnikov - 
de Haas oscillations for the surface electrons in the Bi2X3 family of three-dimensional topological 
insulators (X stands for Te or Se). We also discuss a relation of the obtained results with published 
experimental data on the Shubnikov-de Haas oscillations for this family of topological insulators. 
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Three-dimensional topological insulators (TI) attract 
considerable interest as a new state of solids.^ These ma- 
terials are insulating in the bulk and are metallic on their 
surfaces. Recently, a new important class of these ma- 
terials was discovered that consists of 612X3 compounds 
where X stands for Te or Sei^i^ The angle-resolved pho- 
toemission spectroscopy revealed that the energy bands 
of the surface electrons are practically linear in the wave 
vector k in these compounds, forming a massless Dirac 
cone in their surfaces It is necessary to emphasize 
that although a similar Dirac cone exists in graphene, the 
surface electron states in the TI essentially differ from the 
electron states in graphene. The surface states of the TI 
have no inversion symmetry, and the spin-orbit coupling 
is not weak in these crystals. These features of the TI 
lead to a locking of electron spins with k for the surface 
states, giving rise to Dirac fermions without spin degen- 
eracy. 

Recently, the quantum oscillations of conductivity 
in the magnetic field (the Shubnikov-de Haas oscil- 
lations) were observed in (Bii_a;Sba;)2Se3(^ Bi2Te3fi'^ 
Bi2Te2Se^ia Bi1.5Sbo.5Te1.7Se1.3ri and Bi2Se3^ It was 
found that the observed oscillations result from the sur- 
face electrons of these TI. For such a two-dimensional 
metal the part of the conductivity describing these oscil- 
lations in the semiclassical limit, &(Jxxi has the iorvni^— 
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where Ai are positive amplitudes of the harmonics of the 
periodic in 1/H signal, F and Lp are the frequency and 
the phase of the oscillations, the x-y plane coincides with 
the surface of TI, and the magnetic field H — is 
perpendicular to this surface. The frequency F is de- 
termined by the area S{eF) of the closed orbit of an 
electron with the Fermi energy e^:^ in the space of the 
wave vectors k, 2ttF — ficS{ep)/ef^^'^'^ where e is the 
absolute value of the electron charge. The phase (p is 
given by the constant ip = —2ttj^ that enters the 
semiclassical quantization condition for the electron en- 
ergy in the magnetic field, Eq. ([2]). Both F and (p 
can be experimentally determined by Fourier analysis of 



ScTxxi^/H).^'^ It was found in the experiments^' that 
the phase of the oscillation, ip, practically coincides with 
the appropriate phase for the electrons in graphene, i.e., 
93 w 0. On the other hands, this phase measured in 
other experimental investigations^i^iiii generally differs 
both from the phase for the usual two-dimensional elec- 
tron gas {ip = — tt) and from the phase in graphene. As 
will become clear below, the constant 7 (i.e., ip) depends 
on the so-called Berry phase^^ of electrons and on the 
Zeeman coupling of an electron moment with the mag- 
netic field. In Refs.d, 0, and it is the Zeeman couphng 
that was discussed as a possible source of the nonzero ip. 
Taskin and Andc^^ suggested that a deviation of the dis- 
persion relation e(k) for the surface charge carriers from 
an ideal cone can shift their Berry phase $b from its 
usual value tt characteristic of the Dirac electrons and 
hence can change the phase of the oscillations ip, too. In 
this paper, using the semiclassical approach, we calculate 
both the Berry phase and a contribution of the Zeeman 
coupling to the constant 7 for the surface electrons in the 
Bi2X3 family of TI. We also discuss a relation of these 
results with the published experimental data. 

The semiclassical quantization condition for the energy 
levels e„ of an electron in the magnetic field H reads i^^i^° 
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where n is a large integer. The quantization condition ^ 
is obtainable from the one-band Hamiltonian^^ of a semi- 
classical electron in the magnetic field. This Hamiltonian 
H can be schematically represented in the form: 



^ = £o(k)-^(A.[,') 
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where the subscript marks a two-dimensional electron- 
surface band under study, the term 

M^')(k) = [voxf^o]. (4) 

is proportional to the intraband part of the orbital elec- 
tron moment, Vo = (l/?i)(9eo/9k) is the electron veloc- 
ity, 



f^o(k) 



d 

4,o(i")^Wk,o(r)dr, 



(5) 



2 



(2) 

and /Xg is the sum of the z components of the interband 
orbital moment of the electron and of its spin. In Eq. (O 
the integration is carried out over a unit cell of the crystal 
lattice, and Uk.o(r) is the periodic factor in the Bloch 
wave function of the 0-th band, 

■0k,o = exp(ikr)uk^o- 

Before proceeding to TI, let us briefly review possible 
situations that can occur in crystals and that have an 
effect on the form of Eq. If the spin-orbit interac- 
tion in a crystal is negligible, and if the electron states 
under study have the inversion symmetry, the constant 
7 is given by the formula)^ 



1 
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where 4>b is the Berry phase of the electron orbit F, 
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dk = dn^z Vo]/|vo|; dn is the length of an infinitesimal 
element of the orbit F, and is the unit vector parallel to 
H. It is this situation that occurs in graphene. The part 
<I'_B/27r of Eq. ([6]) results from the term /ig^'' in Hamilto- 

nian ([3]) . The term /Xg^^ is identically equal to zero in this 
case (if one neglects the electron spin). As was shown in 
our paper,— the Berry phase $b of an electron is always 
equal to ±7r when it moves around a Dirac point, i.e., 
if the electron orbit F surrounds this point. Otherwise, 
one has $b = 0. Importantly, the Berry phase does not 
depend on the shape and the size of the electron orbit 
since the inversion and time-reversal symmetries lead to 
rotrio = everywhere except the Dirac point. In other 
words, the result <^b — ±tt remains true at any depen- 
dence eo(k) in the vicinity of the orbit, and 7 can take 
on only the universal values or 1/2. 

If the spin-orbit interaction is not weak, but still there 
is an inversion symmetry for the electron states, all these 
states are double degenerate in spin. In this case one 
has 7 — 1/2, but the quantization rule ^ contains an 
additional term associated with the so-called electron g 
factor, (7, i.e., in this case 7 — >■ 1/2 ± grrisr/Am where m 
is the electron mass and to* is its cyclotron mass. One 

part of this g factor is still determined by the Berry phase, 

(2) 

whereas its second part depends on the term /ig which is 
not vanish now. The complete theory of the g factor for 
electrons in metals was presented in our papers. ^'^■^^ It is 
important that if the strength of the spin-orbit interac- 
tion is small, the g factor approaches a limiting form that 
leads to the same electron spectrum described by Eq. ^ 
with 7 = 1/2 or as in the absence of this interaction^^ 
When the inversion symmetry is absent, formulas for 
the constant 7 can be derived by a siinple modification 
of the expressions given in Refs.l23l andl24l. In particular, 
in the case of TI when the spin-orbit interaction is not 
weak, and hence when the surface electron states are not 



degenerate in the electron spin, we obtain 
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Here the first integral coincides with — $B/27r, while the 

second integral can be considered as a result of the Zee- 

(2) 

man coupling of the electron moment /ig with the mag- 
netic field. For ideal Dirac cone, $s is still equal to ±n 
when the electron orbit surrounds the Dirac pointji^ But 
if the electron band £o(k) is not described by a strictly 
linear dependence in the whole energy interval from the 
Dirac point to the Fermi energy ep, this result for $b 
generally fails, and the Berry phase may differ from tt. It 
is this possibility that was discussed in Ref. [l^. 

The functions Oq (k) and /ig^' (k) can be calculated^^i^^ 
using a k-p Hamiltonian of electrons at a point kg located 
near the electron orbit. This calculation is also possi- 
ble with an effective two-band Hamiltonian that contains 
nonlinear terms in k — kg 1^1^ Such the terms are usually 
introduced into the Hamiltonian in order to take into ac- 
count the other bands different from the two considered 
explicitly. In the magnetic field H the effective two-band 
Hamiltonian should also contain a term linear in H. This 
term is required for the accurate calculation of /ig'^' (k) . 

In Refs. |27| and |2^ the effective two-band Hamiltonian 
was found for the surface electrons of 612X3 family in the 
vicinity of the center of their two-dimensional Brillouin 
zone, 
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where 5, v, c, and gsz are some material-dependent pa- 
rameters; fj-B is the Bohr magneton, fc^ = k^. + ky] 

k± = fca; ± iky] 1 is the unit matrix, and cTi are the 
Pauli matrices. At H = this Hamiltonian describes the 
conduction and valence surface bands. 



ecA^)=bk^±E, 
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which touch each other at the point fc = 0, and which we 
denote by the indexes c and v. Note that these bands are 
not strictly linear functions of k. They also contain the 
quadratic and cubic corrections to the Dirac cone. The 
values of the parameters defining Eq. © for the crystals 
of the 612X3 family were estimated in Ref. [H. 

Using Hamiltonian ^ and formulas (A2), (A3), (A7), 
(A8) of our paper we calculate /^(^^(k) and r2(k), e.g., 
for the conduction band c. 



± 



2E{E 
ckx{k^ - 



E 



Eh 



ikDY 



2m 



(11) 



(12) 



where m is the electron mass, and the signs -I- and — cor- 
respond to {ilc)x and {i^c)y: respectively. Note that Hc 
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is proportional to the parameter c defining the hexagonal 
warping term in Hamiltonian ^ . Without this term, the 
electron moment lies in the x-y plane,— and so it has no 
component along the magnetic field. 

Using formulas (fTTj) . (IT2]) and an expression for the 
velocity Vc — |9ec(k)/9k| that is directly obtainable 
from Eqs. pTI)l . one can find 7, calculating the inte- 
grals over the orbit F in Eq. ([5]). This orbit is de- 
fined by the condition £c(k) — ep where the constant 
ep is the Fermi energy. It is significant that the orbit 

is symmetric relative to the transformation —k^, 

(2) 

ky — > ~ky, while the term /ic is antisymmetric to this 
transformation. Thus, the second integral in Eq. ([8]) is 
always equal to zero. Representing 17 as a sum of the 
symmetric Jl* = [fl(k) + k)]/2 and antisymmetric 
rj" = [Jl(k) — r2(— k)]/2 parts, we find that only the an- 
tisymmetric part gives a nonzero contribution to 7. A 
direct calculation with Eq. yields 
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This rjJJ leads to $5 = — tt, and hence one obtains 7 = 
for any electron orbit in the conduction band (the val- 
ues 7 = 1 and 7 = are equivalent in the semiclassical 
limit). A similar calculation gives $b = tt, and 7 = for 
the valence band. Interestingly, if the Dirac cone in the 
Brillouin zone were shifted from the point k = or were 
deformed asymmetrically, the Berry phase and the Zee- 
man term in Eq. (|8]) would not have the universal values 
and would depend on e^; see below. 

In fact, the obtained result, 7 = 0, is exclusively caused 
by the time-reversal symmetry of the Hamiltonian. This 
symmetry alone dictates that any electron orbit is sym- 
metric relative to the transformation k — > — k, and that 
the Hamiltonian should have the form 



i=l 
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+ ^i/,[5o(k)i+5^5.(k)<T.], (14) 

where arbitrary functions /io(k), /ii(k) and .go(k), .9i(k) 
satisfy only the requirements: /io(k) = /io(— k), 5o(~k) = 
-go(k), and h,{k) = -h,{-k), gt{-k) g,{k). With 
Eq. we find the following expressions for r2c(k) and 

He (k) generalizing formulas ([TT|) - p^ : 
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where Siji is the completely antisymmetric unit tensor 
with £123 = 1. A direct calculation shows that rotfJ" = 
everywhere except the Dirac point for any functions 
/ii(k), /i2(k). Hence, for the calculation of the Berry 
phase, it is sufficient to consider Jp fJ^dk in the imme- 
diate vicinity of the Dirac point where these functions 

reduce to linear ones. For linear functions /ii(k), /i2(k) 

{2.) 

one indeed obtains $b = ±7r. As to the term /ic (k), it 
is antisymmetric relative to the transformation k — > — k 
and does not contribute to 7. However, we emphasize 
here that formulas ([H]) - ([TT]) are, in principle, applica- 
ble to a general case (i.e., not only to the Bi2X3 family of 
TI) , and that these formulas can lead to nontrivial values 
of 7 if the appropriate electron orbits in some TI becomes 
asymmetric with respect to the point k = 0. 

The asymmetry of the orbit can appear when the ex- 
ternal magnetic field is not strictly perpendicular to the 
surface of TI, i.e. when it has an in-plain component H|| . 
In this case, e.g., Hamiltonian ^ contains an additional 
terro^l 
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where is the in-plain g factor, (jy ~ li^ This term 
shifts the Dirac cone of the surface electrons in their 
Brillouin zone from the point k = to the point 
k* — 5|| (MB/2'y)[H X iz],^"^ and the electron orbits be- 
come asymmetric. This asymmetry generally leads to 
a nonzero value of 7. However, in the first order in the 
small parameter g^^iiBH^^/ep the correction to 7 vanishes 
for the case of Hamiltonian ([9]) since the angular depen- 
dence cos3(/) characteristic of the factor ft.3(k) in Eq. 
differs from the angular dependence cos (f) (or sin 0) that 
describes the deformation of the electron orbits due to the 
field Here 4> is the angle between k and the x axis. 
Due to this difference in the angular dependences, the 
appropriate integrals over (j) determining the correction 
vanish. Thus, for the family Bi2X3 of TI the effect of the 
in-plain magnetic field on 7 is negligible. This conclusion 
is in an agreement with the experimental data-' show- 
ing that the phase of the Shubnikov - de Haas oscillations 
does not alter after applying the field ffy. 

We now discuss the relation of the obtained results 
with the published experimental data. In Ref. ^ it was 
shown that in (Bii_a;Sba;)2Se3 the electron Landau levels 
e„ at small n can be well described by the formula 
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where gsz ~ 50. This formula gives the exact Landau 
levels of Hamiltonian ^ a.t b = c — 0. Let us rewrite 
expression (jl9l) in the form of quantization condition ^ , 
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where S'(e) = Trfc^ = TT{e/v)'^ is the area of the electron 
orbit in the k-space, and consider the case of large n 
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(i.e., of relatively small H). This case just corresponds to 
the seniiclassical approximation used in our calculations. 
Now the second term in the right hand side of formula 
(PU)) is small as compared to the first one, and we may 
express H via n as follows: H « hcS/ (2TTen). Inserting 
this expression for H in the second term of Eq. (|20|). we 
arrive at 
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Strictly speaking, the right hand side of the quantiza- 
tion condition is a series in \/n [i.e., it has the form 
n + const + (1/n.) + ...] in which only the first two terms, 
71 + 7, are usually kept. It is this cut of the series that 
yields Eq. Formula (PT|) shows that spectrum 
in the semiclassical limit (n ^ 1) leads to 7 = in an 
agreement with our results obtained above and that in 
contrast with the case of the usual electron g finite 
Qsz in Eq. (jl9p gives a contribution only to the term of 
the order of \/n. In other words, the nonzero values of 
7 found in the experiments,—'^'^ are due to the difference 
of the exact spectrum from the semiclassial one at not 
too large values of n. This conclusion is supported by an 



analysis of the Shubnikov - de Haas oscillations carried 
out by Taskin and AndoiS. for various TI. Note that in 
graphene one has a relatively small value of gsz ^ 2, and 
the semiclassical spectrum defined by Eq. © practically 
coincides with the exact spectrum (fT9|) even at = 
and 1. The large value of gsz ~ 50 for (Bii_2;Sba;)2Se3 is 
probably due to a relatively small gap between the Dirac- 
point energy and the energies of some other surface bands 
at k = for this TI. 



In summary, within the semiclassical approach, we de- 
rive formula for the constant 7 that enters the quan- 
tization condition ([2|) and defines the phase — — 27r7 
of the Shubnikov - de Haas oscillations in TI. This 7 is 
determined by the Berry phase of an electron orbit and 
by a part of the electron moment averaged over the or- 
bit. Since the Dirac point lies at k = for 612X3 family 
of TI, we find that the Berry phase is equal to tt, while 
the averaged moment is zero, and hence one always has 
7 = for this family of TI. The nonzero values of 7 found 
in some experiments appear to be due to the difference 
of the exact electron spectrum from the semiclassial one 
at not too large quantizing numbers n. 
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